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ABSTRACT 

cn 
o . 

O ■ We propose a model to explain how a Gamma Rays Burst can take place days 

^ or years after a supernova explosion. Our model is based on the conversion of a 

D ■ pure hadronic star (neutron star) into a star made at least in part of deconfined 

quark matter. The conversion process can be delayed if the surface tension at 

the interface between hadronic and deconfined-quark-matter phases is taken into 

account. The nucleation time {i.e. the time to form a critical-size drop of quark 

^ ■ matter) can be extremely long if the mass of the star is small. Via mass accretion 

^ ■ the nucleation time can be dramaticaly reduced and the star is finally converted 

CNj . into the stable configuration. A huge amount of energy, of the order of 10^^- 

0\' 1 n53 



o 



10 erg, is released during the conversion process and can produce a powerful 



CN . Gamma Ray Burst. The delay between the supernova explosion generating the 

. metastable neutron star and the new collapse can explain the delay proposed in 

GRB990705 (Amati et al. 2000) and in GRB011211 (Reeves et al. 2002). 



O 

! of state 
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1. Introduction 

The discovery of a transient (13 s) absorption feature in the prompt emission of the 
40 s Gamma Ray Burst (GRB) of July 5, 1999 (GRB990705) (Amati et al. 2000) and 
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the evidence of emission features in the afterglow of several GRBs (Piro et al. 1999; Yoshida 
et al. 1999; Piro et al. 2000; Antonelli et al. 2000; Reeves et al. 2002) have stimulated the 
interpretation of these characteristics in the context of the fireball model of GRBs. Amati et 
al. (2000) attribute the transient absorption feature of GRB990705 (energy released ~ 10^^ 
erg assuming isotropy) to a redshifted K edge of Iron contained in an environment not far 
from the GRB site (~ 0.1 pc) and crossed by the GRB emission. They estimate an Iron 
abundance typical of a supernova (SN) environment (Ape ~ 75) and a time delay of about 
10 years between the SN explosion and the GRB event. Lazzati et al. (2001) give a different 
interpretation of the absorption feature, in terms of a redshifted resonance scattering feature 
of H-hke Iron (transition ls-2p, E'rest = 6.927 keV) in an inhomogeneous high-velocity 
outflow, but invoke a Iron rich environment as well, due to a preceding SN explosion, even 
if a shorter time delay (~ 1 yr) between SN and GRB is inferred. A SN explosion preceding 
the GRB event is also inferred for explaining the properties of the emission features in the 
X-ray afterglow spectrum of GRB000214 (Antonelh et al. 2000) and GRB991216 (Piro et 
al. 2000). In the latter case it cannot be excluded that the SN explosion occured days 
or weeks before the GRB (Rees & Meszaros 2000). Reeves et al. (2002), to explain the 
multiple emission features observed in the afterglow spectrum of GRB011211 (time duration 
of ~ 270 s, isotropic gamma-ray energy of 5 x 10^^ erg), invoke a SN explosion preceding 
the GRB event by ~ 4 days (in the isotropic limit, a minimum of 10 hrs). Even if other 
interpretations for the afterglow emission lines are possible which do not involve a previous 
SN explosion (e.g., (Rees & Meszaros 2000; Meszaros & Rees 2001)), this explosion seems to 
be the most likely way to explain the transient absorption line observed from GRB990705 
(Bottcher et al. 2002). In conclusion, the previous observations suggest that, at least for a 
certain number of GRBs, a SN explosion happened before the GRB, with a time interval 
between the two events ranging from a few hours to a few years. In this context, an attractive 
scenario is that described by the supmnova model (Vietri & Stella 1998) for GRBs. In this 
model, the GRB is the result of the collapse to a black hole (BH) of a supramassive fast 
rotating neutron star (NS), as it loses angular momentum. According to this model the NS 
is produced in the SN explosion preceding the GRB event. The initial barionic mass Mb of 
the NS is assumed to be above the maximum baryonic mass for non-rotating configurations. 
However, as also noticed by Bottcher et al. (2002), on the basis of realistic calculations of 
collapsing NS (Fryer & Woosley 1998), in these collapses too much baryonic material is 
ejected and thus the energy output is expected to be too small to produce GRBs. Even if 
the introduction of magnetic fields or beaming could overcome this limitation, in any case, 
the GRB duration from a NS collapse should be very short (^ 1 s), much shorter than that 
observed from GRB990705. 

In this Letter, we propose an alternative model to explain the existence of GRBs as- 
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sociated with previous SN explosions. In this model, unlike the supranova model, the NS 
collapse to BH is replaced by the conversion from a mctastable, purely hadronic star (neu- 
tron star) into a more compact star in which deconfined quark matter (QM) is present. This 
possibility has already been discussed in the literature (Cheng & Dai 1996; Bombaci & Datta 
2000; Wang et al. 2000; Ouyed ct al. 2002). The new and crucial idea we introduce here, is 
the metastability of the purely hadronic star due to the existence of a non-vanishing surface 
tension at the interface separating hadronic matter from quark matter. The mean-life time 
of the metastable NS can then be connected to the delay between the supernova explosion 
and the GRB. As we shall see, in our model we can easily obtain a burst lasting tens seconds, 
in agreement with the observations. The order of magnitude of the energy released is also 
the appropriate one. 



2. Quark Matter nucleation in compact stars 

Recently various possibilities have been discussed in the literature to get compact stars 
in which matter is, partially or totally, in a state of deconfined quarks (see e.g. Glendenning 
(2000); Heiselberg & Hjorth- Jensen (2000); Drago & Lavagno (2001)). Concerning the stellar 
quark content, it is possible to have three different classes of compact stars: a) purely 
hadronic stars (HS), in which no fraction of QM is present; b) hybrid stars (HyS), in which 
only at the center of the star QM is present either as a mixed phase of deconfined quarks 
and hadrons or as a pure phase; c)quark stars (QS), in which the surface of the star is made 
of matter having a large density, of the order of nuclear matter saturation density or larger, 
and the bulk of the star is made of deconfined QM. The sizeable amount of observational 
data collected by the new generations of X-ray sateUites, has provided a growing body of 
evidence for the existence of very compact stars, which could be HyS or QS (Bombaci 1997; 
Cheng et al. 1998; Li et al. 1999a,b; Xu 2002; Drake et al. 2002). 

In our scenario, we consider a purely HS whose central density (pressure) is increasing 
due to spin-down or due to mass accretion (e.g., from fallback of ejected material in the SN 
explosion). As the central density approaches the deconfinement critical density, a virtual 
drop of quark matter can be formed in the central region of the star. The fliictuations of a 
spherical droplet of quark matter having a radius R are regulated by a potential energy of 
the form (Lifshitz & Kagan 1972) 

4 

U{R) = -TTR\{i^q - Hh) + 4:TTaR^ + 87r-fR (1) 
3 

where Ug is the quark baryon density, /ih and /ig are the hadronic and quark chemical 
potentials at a fixed pressure P, and a is the surface tension for the surface separating 
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quarks from hadrons. Finally, the term containing 7 is the so called curvature energy. The 
value of the surface tension a is poorly known, and typical values used in the literature range 
from 10 to 50 MeV/fm^ (Heiselbcrg ct al. 1993; lida & Sato 1998). Following the work of lida 
&; Sato (1998), we have assumed that the term with a takes into account in an effective way 
also the curvature energy. The term with 7 is discussed e.g. by Masden (1993), while other 
more complicated terms, connected with the Coulomb energy, are discussed in the literature 
(Heiselberg et al. 1993; lida & Sato 1998). We have neglected them in our analysis since 
they do not dramatically modify both the nucleation time and the energy associated with 
the transition into the stable quark matter configuration. 

If the temperature is low enough, the process of formation of a bubble having a critical 
radius proceeds through quantum tunnelling and it can be computed using a semiclassical 
approximation. The procedure is rather straightforward. First one computes, using the semi- 
classical (WKB) approximation, the ground state energy Eq and the oscillation frequency 
vq of the virtual QM drop in the potential well U{R). Then it is possible to calculate in a 
relativistic frame the probability of tunneling as (lida & Sato 1998): 

po = exp[-^] (2) 

where ^ 

A{E) = 2 [ ^ dR^J[2M{R) + E- U{R)][U{R) - E] . (3) 

JR- 

Here R± are the classical turning points and 

M(i?) = 47rp,(l - ^)2i?3 , (4) 

Ph being the hadronic energy density (here and in the following we adopt the so-called 
"natural units", in which h — c — 1). Uh and Ug are the baryonic densities at a same and 
given pressure in the hadronic and quark phase, respectively. The nucleation time is then 
equal to 

T = (uoPoN,)-^ , (5) 

where is the number of centers of droplet formation in the star, and it is of the order of 
10^8 (lida & Sato 1998). 



3. Results 



The typical mass-radius relations for the three types of stars we are discussing can be 
found e.g. in Fig. 3 of Drago & Lavagno (2001), where a relativistic non-hnear Walecka-type 
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model (Glendenning & Moszkowski 1991) has been used to describe the hadronic phase. 
As it appears, stars containing QM (either HyS or QS) are more compact than purely HS. 
In particular QS can have much smaller radii then HS when they have a small mass. In 
our scenario a metastable HS having a mass of, e.g., 1.3 Mq and a radius of ~ 13 km can 
collapse into an HyS having a radius of 10.5 km or into a QS with radius ~ 9 km (with 
respect to the results of Fig. 3 of Drago & Lavagno (2001), here we use an equation of state 
(EOS) which includes hyperonic degrees of freedom). The nature of the stable configuration 
reached after the stellar conversion (i.e. an HyS or a QS) will depend on the parameters of 
the quark phase EOS. 

The time needed to form a critical droplet of deconfined quark matter can be calculated 
for different values of the stellar central pressure Pc (which enters in the expression of the 
energy barrier in eq. (1)) and it can be plotted as a function of the gravitational mass Mhs of 
the HS corresponding to that given value of central pressure. The results of our calculations 
for a specific EOS of hadronic matter (the GM3 model with hyperons of Glendenning & 
Moszkowski (1991)) are reported in Fig. 1, where each curve refers to a different value of 
the bag constant B. If we assume, for example, B^^'^ = 170 McV (which corresponds to 
B = 109 MeV/fm^) and the initial mass of the HS to be Mhs = 1-32 Mq, we find that the 
"life time" for this star is about 10^^ years. As the star accretes a small amount of matter, 
the consequential increase of the central pressure lead to a huge reduction of the nucleation 
time, and, as a result, to a dramatic reduction of the HS life time. For our HS with initial 
mass of 1.32 Mq the accretion of about 0.01 Mq reduces the star hfe time to a few years. 
We would hke to stress that in our model the delay between the SN explosion and the GRB 
is regulated by the mass accretion rate, rather then by the mass and the spinning of the 
metastable star itself. Since the mass accretion rate is generally larger during the first days 
after the SN explosion, a delay of a few days will be rather typical in our scenario. However, 
longer delays are also possible if the material ejected during the SN explosion has a small 
fallback. 

In the model we are presenting, the GRB is due to the cooling of the justly formed 
HyS or QS via neutrino-antineutrino emission (and maybe also via emission of axion-like 
particles, see below). The subsequent neutrino-antineutrino annihilation generates the GRB. 
In our scenario the duration of the prompt emission of the GRB is therefore regulated by 
two mechanisms: 1) the time needed for the conversion of the HS into a HyS or QS, once 
a critical-size droplet is formed and 2) the coohng time of the justly formed HyS or QS. 
Concerning the time needed for the conversion into QM of at least a fraction of the star, 
the seminal work by Olinto (1987) has been reconsidered by Horvath & Benvenuto (1988). 
The conclusion of this latter work is that the stellar conversion is a very fast process, having 
a duration much shorter than Is. On the other hand, the neutrino trapping time, which 
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provides the cooling time of a compact object, is of the order of a few ten seconds (Prakash 
et al. 1997), and it gives the typical duration of the GRB in our model. In Table 1 we 
give the measured duration and the estimated electromagnetic energy (assuming isotropic 
emission) of the GRBs associated with Fe emission or absorption lines. All bursts last at 
least 10 s. According to our model the firsts few ten seconds correspond to a prompt 7- 
rays emission, while the subsequent emission should be interpreted as the beginning of the 
afterglow. Actually it has been found that at least the second half of the prompt emission 
of long bursts is likely due to afterglow (Prontera et al. 2000). We would like to remark 
however that we are not suggesting that all the GRBs should be explained in our model. In 
particular long and energetic bursts could be originated e.g. by coUapsars (MacFadyen & 
Woosley 1999). On the other hand, the variety of GRB durations could be explained within 
the QS formation scenario itself making use of the "unstable photon decay" mechanism 
proposed by Guyed & Sannino (2002). 

Next we consider the total energy A£' released in the transition from a metastable 
HS (with hyperonic degrees of freedom) to HyS or QS (which final state is reached in this 
transition depends on the details of the QM EOS and in particular on the value of the bag 
constant). The energy released is calculated as the difference between the gravitational mass 
of the metastable HS and that of the final stable HyS (or QS) having the same baryonic 
mass (Bombaci & Datta 2000). In Table 2 we report the energy released for various values 
of the bag constant B and of the surface tension a. Notice that the transition will take place 
when the nucleation time will be reduced to a value of the order years, due e.g. to mass 
accretion on the HS (recall the exponential dependence of the nucleation time on the mass 
of the HS, as shown in Fig. 1). Therefore the total energy released in the collapse will be 
always of the same order of magnitude, once the parameters of the model have been fixed. 
As shown in Table 2, the released energy is in the range (3 — 5) x 10^^ ergs for all the sensible 
choices of the EOS parameters. The "critical mass" Ma of a metastable HS having a lifetime 
T = 1 yr is in the range (0.9 — 1.4) Mq. When the mass of the HS reaches a value near 
Mcr, the conversion process takes place. It is worth mentioning that the energy released in 
the conversion can be larger if a diquark condensate forms inside QM, see e.g. Hong et al. 
(2001). 

To generate a strong GRB, an efficient mechanism to transfer the energy released in the 
collapse into an electron-photon plasma is needed. In an earlier work (Fryer & Woosley 1998) 
it was this difficulty that hampered the possibility to connect GRBs and the hadronic-quark 
matter phase transition in compact stars. Only more recently it was noticed (Salmonson 
& Wilson 1999) that near the surface of a compact stellar object, due to general relativity 
effects, the efficiency of the neutrino-antineutrino annihilation into e~^e~ pairs is strongly 
enhanced with respect to the Newtonion case. The efficiency of the conversion of neutrinos 
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in e^e~ pairs could be as high as 10%. In the computation of the energy associated with 
the final GRB we must take into account the possibility of a moderate anisotropy of the 
electron motion, due to the presence of the magnetic field of the star, which will in turn 
generate a moderate anisotropy of the burst emission ^. Other anisotropics in the GRB 
emission could be generated by the rotation of the star which could affect the efficiency 
of the neutrino-antineutrino annihilation due to general relativity effects. On the basis of 
these considerations, the energy deposited in the burst could be sufficient to explain the 
isotropic energy of the GRBs fisted in Table 1. We must also recall that more efficient ways 
to generate photons and/or e^e" pairs have been proposed in the literature, based on the 
decay of axion-like particles (Berezhiani & Drago 2000). This mechanism would have an 
extremely high efficiency and would transfer most of the energy produced in the collapse 
into GRB electromagnetic energy. 

There are various specific signatures of the mechanism we are suggesting. First, two 
classes of stars having similar masses but rather different radii should exist: a) pure (metasta- 
ble) HS, with radii in the range 12-20 km, as is the case of the compact star IE 1207.4-5209, 
assuming M = IAMq (Sanwal et al. 2002), and b) HyS or QS with radii in the range 6-8 
km (Bombaci 1997; Li et al. 1999a,b; Drake et al. 2002). Second, all the GRBs generated by 
the present mechanism should have approximately the same isotropic energy and a duration 
of at least 10 s. 



4. Conclusions 

We propose the following origin for at least some of the GRBs having a duration of 
tens of seconds. They can be associated with the transition from a metastable HS to a 
more compact HyS or a QS. The time delay between the supernova explosion originating the 
metastable HS and the GRB is regulated by the process of matter accretion on the HS. While 
most of the stellar objects obtained by a SN explosion will possibly have a mass larger than 
Mcr and will therefore directly stabilize as HyS or QS at the moment of the SN explosion, in 
a few cases the mass of the protoneutron star will be low enough not to allow the immediate 
production of QM inside the star. Only when the star will acquire cnoiigh mass, the process 
of QM formation could take place. Due to the surface tension between the hadronic matter 
and the QM the star will become metastable. The later collapse into a stable HyS or QS will 
generate a powerfuU GRB. It can be interesting to notice that, in order to have a not too 



^Dramatic effects of a time-dependent magnetic field have been discussed e.g. by Kluzniak & Ruderman 
(1998). 
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small value for Mcr, a relatively large value for the bag constant B has to be choosen, S^/^ ~ 
170 MeV, which turns out to be the prefered value in many hadronic physics calculations 
(see e.g. (Steffens et al. 1995)). In this situation the final state is an HyS and not a QS. 

It is a pleasure to thank Elena Plan and Luciano RezzoUa for very useful discussions. 
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Fig. 1. — Time needed to form a quark matter droplet as a function of the mass of the HS 
for five different values of S^/^[MeV]. The hadronic phase has the GM3 parameters set with 
hyperons, the quark phase has = 150 MeV and a surface tension a — 30 MeV/fm^ is 
assumed. 
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Table 1: Duration and energy released (assuming isotropy) of the GRB associated with the 
presence of emission or absorption Fe lines in the spectrum. The data have been estracted 
from Amati et al. (2002) and Bloom et al. (2001). 
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BV4[MeV] 


(7[MeV/fm2] 




A E [10^^ erg] 


170 


20 


1.25 


30.0 


170 


30 


1.33 


33.5 


170 


40 


1.39 


38.0 
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30 


1.15 


38.6 


160 


30 


0.91 


45.7 



Tabic 2: Critical mass Mcr of the metastable hadronic star (in unit of the mass of the sun 
Mq = 1.989 X lO-'^-^ g) and energy released AE in the conversion to hybrid star assuming 
the hadronic star mean life time r equal to 1 year. Results are reported for various choices 
of the surface tension a and of the bag constant B. The strange quark mass is taken equal 
to 150 MeV. For the hadronic matter EOS the GM3 model with hyperons (Glendenning & 
Moszkowski 1991) has been used. 



